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Abstract 

In the present paper the transport of impurities driven by trapped electron (TE) mode 
turbulence is studied. Non-linear (NL) gyrokinetic simulations using the code GENE are 
compared with results from quasilinear (QL) gyrokinetic simulations and a computationally 
efficient fluid model. The main focus is on model comparisons for electron temperature gra- 
dient driven turbulence regarding the sign of the convective impurity velocity (pinch) and 
the impurity density gradient R/ Ln^ (peaking factor) for zero impurity flux. In particular, 

P^ the scaling of the impurity peaking factors with impurity charge Z and with driving temper- 

ature gradient is investigated and compared with the results for Ion Temperature Gradient 
(ITG) driven turbulence. In addition, the impurity peaking is compared to the main ion 

I peaking obtained by a self-consistent fluid calculation of the density gradients corresponding 

^— ( to zero particle fluxes. 

^^ For the scaling of the peaking factor with impurity charge Z, a weak dependence is 

fvg obtained from NL GENE and fluid simulations. The QL GENE results show a stronger 

dependence for low Z impurities and overestimates the peaking factor by up to a factor of 

two in this region. As in the case of ITG dominated turbulence, the peaking factors saturate 

as Z increases, at a level much below neoclassical predictions. However, the scaling with Z 

is weak or reversed as compared to the ITG case. 

ri The scaling of impurity peaking with the background temperature gradients is found to 

ri be weak in the NL GENE and fluid simulations. The QL results are also here found to 

j^ significantly overestimate the peaking factor for low Z values. 
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For the parameters considered, the background density gradient for zero particle fiux is 
found to be slightly larger than the corresponding impurity zero fiux gradient. 



The transport properties of impurities is of high relevance for the performance and optimisation 

of magnetic fusion devices. For instance, the possible accumulation of He ash in the core of the 

reactor plasma will serve to dilute the fuel, thus lowering fusion power. Heavier impurity species, 

^kJ originating from the plasma-facing surfaces, may also accumulate in the core, and wall-impurities 

/«~s of relatively low density may lead to unacceptable energy losses in the form of radiation. In an 

00 operational power plant, both impurities of low and high charge numbers will be present. 

OO In the confinement zone of tokamaks, the transport of the background species is usually 

^^ dominated by turbulence. The Trapped Electron (TE) mode and the Ion Temperature Gradi- 

f~^ ent (ITGpJ mode are expected to be the main contributors. Turbulent impurity transport has 

^D been investigated in a number of theoreticall^tSl] g^j^^^ experimentall^^^ papers. In tokamak ex- 

] periments, also the impurity transport is usually dominated by turbulence, resulting in impurity 

. . peaking factors well below the neoclassical predictions, i''' 22 27 28 rp-|_^g main theoretical effort has, 

^ with a few exceptions,! ^ * ^ * ^" * ^ hitherto been devoted to quasilinear studies, primarily focused on 

ITG mode driven impurity transport. For the directly reactor relevant regimes, however, where 

; ^ a-particle heating dominates, as will be the case in the ITER device, or in electron cyclotron 

Cd resonance heated plasmas, TE mode driven impurity transport will likely be important. 



^^also commonly referred to as the r]i mode 



In the present study, transport of impurities driven by TE mode turbulence is investigated by 
NL gyrokinetic simulations using the code GENE.I3SH31I xhe simulation results are compared to 
QL gyrokinetic simulations as well as results obtained from a multi- fluid model.'^The fluid model 
is employed for the dual purposes of benchmarking a computationally efficient model, suitable 
for predictive simulations, and interpreting the results. The TE mode results are compared with 
the more well known results for ITG mode dominated turbulence, obtained from QL gyrokinetic 
and fluid simulations. 

The impurity diffusivity (Dz) and convective velocity (Vz) are estimated from simulation 
data, and from these the zero-flux peaking factor impurity density gradient {R/ L^^ = —RVz/Dz), 
also referred to as the impurity peaking factor (PF) is derived. This quantity expresses the im- 
purity density gradient at which the convective and diffusive transport of impurities are exactly 
balanced. The sign of PF is of special interest, as it determines whether the impurities are 
subject to an inward {PF > 0) or outward {PF < 0) pinch. Scalings of peaking factors with im- 
purity charge {Z), electron and ion temperature gradients (VTe^j), and electron density gradient 
(Vrig), are studied, giving particular attention to VTg driven TE mode impurity transport. The 
results are compared and contrasted with results from previous studies focused on ITG driven 
impurity transport. In addition, the impurity peaking relative to the main ion peaking in the 
plasma core obtained from a self-consistent treatment of the particle fluxes will be discussed. 
This section is experimentally relevant in situations with edge particle fuelling where the steady 
state gradient corresponds to zero particle flux. 

The remainder of the paper is structured as follows: first the transport models are reviewed. 



beginning with of the fluid model employed (section II. A I where the focus is on the impurity 



dynamics. This is followed by a brief introduction of the gyrokinetic model and the GENE code 



(section II. B I and a section on the specifics of the simulations (section III). After this, the main 



results are covered, including discussion and interpretation of the acquired results (section IV I. 
The final section of the paper is a summary of the main conclusions to be drawn (section |V[). 

II Transport models 
II. A Fluid theory 

The Wciland multi-fluid modell^Sl consists of coupled sets of equations for each constituent particle 
species: main ions, electrons and impurities.l ^ l ^^ l -'^^ l ^^ l '^'^ l Effects of toroidal rotation is not included 
here. Neglecting finite Larmor radius effects, the impurity equations of continuity, parallel motion 
and energy take the form of Eq. ([T])-(|3|. 

{uJ + Tz)nz- \— \](t) + TzTz ==0 (1) 



k\\5v\\7 Z ~ Tv / ~ \ 

In Eq. (IT])-(l3|, nz — Suz/nz is the density, (j) = ecfy/Te the electrostatic potential, Tz ~ 5Tz/Tz 
the temperature and 5v\\z the parallel velocity. The normalised eigenvalue and wave vector of 
the eigenmodes are a; = w^ + 17 and k = k\\z + k±, ~denoting normalisation with respect to 
the electron magnetic drift frequency uJoe- The normalised scale lengths can be assumed to 
be constant for the flux tube domain considered, and are defined as -r^ = —-S^-a^. where 

' Lx ^i or ^ 



R is the major radius of the tokamak, Xj — rij, Tj for species j. The other parameters are 
defined as follows: r| = XTz/ZTf. with A = cos9 + sOsinO for the poloidal angle 9, tz = Tz/Te, 
Az — mz /mi « 2Z is the impurity mass number, Z is the impurity charge. Further, s is the 
magnetic shear and q* = 2qkgps, where q is the safety factor, ps = Cg/fict is the ion sound scale 
with the ion sound speed Cs = yjT^jrai and the ion cyclotron frequency fid — eB/rrii. Effects 
of curvature enter the equations through the magnetic drift, defined as ujdz = ^rT'^^i^)' which 
originates from the compression of the E x B drift velocity, the diamagnetic drift velocity and 
the diamagnetic heat fiow. Curvature effects from the stress tensor enter as 2r^ at the left hand 
side of Eq. ^. 

Combining Eq. ([T])-(|3]), while neglecting pressure perturbations in Eq. (pi for simplicity, the 
relation of the electrostatic potential and impurity density nz becomes: 
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The main ion and electron response is calculated from the corresponding fluid equations for 
ions and trapped electrons. The electron response is given by a trapped and a free part such 
that ^^ = ft^^ + (1 — /t) 5^1 i-e. the free electrons are assumed to be adiabatic and thus to 
follow the Boltzmann distribution: 611^^/11^^ = e^/Tg. 

The equations are closed by the assumption of quasi-ncutrality: 

= (l-Z/z) \-Zfz , (6) 

He 111 nz 

where /z = ^^ is the fraction of impurities. 

Thus an eigenvalue equation for TE and ITG modes is obtained in the presence of impurities. 
Assuming a strongly ballooning eigenfunction withlSS k?: = (Sq^R^) , the eigenvalue equation 
is reduced to a system of algebraic equations that is solved numerically. The sensitivity of the 
fluid results ti the choice of fcn will be examined in section |IV. Al below. 



The zero-flux impurity peaking factor, defined as PF ~ ^^o for the value of the impurity 
density gradient that give zero impurity flux, quantifles the balance of convective and diffusive 
impurity transport. Its derivation relics on the fact that the transport of a trace impurity species 
can be described locally by a diffusive and a convective part. In the trace impurity limit, i.e. for 
Z fz — )■ in Eq. [6] the impurity flux Fz becomes a linear function of Vn^, offset by a convective 
velocity or "pinch" Vz- The resulting expression can be seen in Eq. ([7|, where nz is the density of 
the impurity species and R is the major radius of the tokamak, and both the diffusion coefficient 
{Dz) and the convective velocity {Vz) are independent of Vnzf^ Setting /^ = in Eq. (It]) 
yields the interpretation of PF as the gradient of zero-impurity flux. 

Fnz=-Dz^nz + nzVz^ — ^Dz^ + RVz, (7) 

nz Fnz 

where the right hand side of the equivalence is arrived at by assuming —Vuz = 1/ Lnz- The 
relationship of PF to Dz and Vz is illustrated in Fig. [T] 

The impurity particle flux at the left hand side of Eq. (l7| can be written as: 



r„2 = (SuzVexb) = -nzPsCs ( nz 
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(8) 



The angled brackets imply a time and space average over all unstable modes. Performing this 
averaging for a fixed length scale kgps of the turbulence, the following expression is reached: 
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where Ni =uj — 2t^ is introduced. 

In the following it is assumed that the turbulence is isotropic in the radial and poloidal 
directions (r and respectively; krPs = kgps), with a saturated fluctuation level of 1 0^1 = 
--^jr-r — .IS^l A brief review of the different mechanisms responsible for the impurity transport, 
as identified in previous studiesHEE3 is given here. The first term in Eq. (l9| corresponds to the 
diffusive part of Eq. (l7| , whereas the three subsequent terms correspond to the convective part 
of the transport of the impurity species. Of this, the R/Ltz term is the thermodiffusion, the 
sign of which is governed mainly by the real frequency, w^. For TE modes, cD,. > 0, and for 
ITG modes w^ < 0, resulting the thermodiffusion generally giving an inward pinch for TE modes 
and an outward pinch for ITG modes. Due to the Z-dependence in r^, this term scales as 
Vz"^ ^ {^1 Z){R/ Ltz) to leading order, rendering it unimportant for large Z impurity species, 
but it is important for lighter elements, such as the Helium ash. Further, the (A) term gives 
the curvature pinch, which is usually inward, and the final term is the parallel compression term 
for the impurities. As opposed to the thermodiffusion, the parallel compression pinch is usually 
outward for TE modes and inward for ITG modes. Its Z dependence is V^ ^ Z/Azk?^ ^ Z/Azq^, 
but since Az ~ 2Z this is is expected to be a very weak scaling.E^ Effects of toroidal rotation on 
the impurity transport has recently been studied ,^^1211 but will not be considered here. 



II.B Gyrokinetics - the GENE code 

The GENE code'22tl22136| jg a massively parallel gyrokinetic Vlasov code, solving the nonlinear 
time evolution of the gyrokinetic distribution functions on a fixed grid in phase space. The 
gyrokinetic equations are derived from the kinetic equations by performing an average of the 
particles' gyrations around the field lines, so that the equations follow the centre of gyration, 
rather than the explicit orbits .I^SEnixhis reduces the velocity space coordinates from three to two 
directions: parallel and perpendicular velocity. Following the conventions of GENE, these are 
referred to as v and /i respectively. In real space, the radial (x) and bi-normal (y) dependencies 
are treated spectrally, i.e. those directions are discretised explicitly in fc-space, whereas the 
toroidal (z) direction is discretised in real space. Because all phase space coordinates are coupled 
nonlinearly, the decrease from six to five phase space coordinates means a significant increase 
in computational efficiency. This simplification of the equations is appropriate if the gyro-radii 
are small compared to the turbulent features, and the cyclotron frequencies large compared to 



the frequencies of the turbulent phenomena, conditions which generally hold in core tokaniak 
plasmas.l^2Hini 

In this paper, GENE simulations are performed in a flux tube geometry with periodic bound- 
ary conditions in the perpendicular directions. The flux tube is in essence a box that is elongated 
and twisted along with the B field as the field lines traverse the tokamak. Its application relies 
on the assumption that the scales of the phenomena of interest are all small compared to the size 
of the whole plasma. This is generally held to be true in the core of the plasma. A cross-section 
of the fiux tube is presented in Fig. |3] There the size of the turbulent features can be seen, and 
a comparison of their size to the flux tube's perpendicular resolution of ~ 125 x 125 main ion 
gyro-radii and the overall box size indicates that the resolution and flux tube dimensions are 
adequate; see Section |TTT] for more details on how the resolution was chosen. 

The data presented in Fig. [3] is computed from the raw field data. By integrating further, 
scalar quantities can be obtained; those quantities are often the most interesting from a physics 
perspective, since they are easier to compare both to theoretical, experimental, and other numer- 
ical results. In this study, the scalar impurity fiux Fz is of most interest. Time series showing 
the fiuctuations in the main ion density and the impurity fiux for a nonlinear GENE simulation 
are presented in Fig. [2] 

GENE can also be run in quasilinear mode, a method that is considerably less demanding 
when it comes computer resources, since the non- linear coupling between length scales is ig- 
nored.^ '^^ l '^^ l The method used here only captures the contribution from the most unstable, not 
sub-dominant modes, and only for the particular length scale kgps of choice. If the length scale 
is chosen appropriately, however, the quasilinear simulation will capture the essential features of 
the transport mechanism, and it is useful for getting a qualitative understanding of the physical 
processes. 

Ill Simulations 

In this paper, the transport of impurities has been studied numerically, by calculating the im- 
purity peaking factor {PF) for impurities with various impurity charge (Z) and varying values 
of the driving background gradients. The process of calculating the peaking factor is illustrated 
in Fig. IT] The impurity particle fiux Fz is computed for Vn^ in the vicinity of Fz = 0, taking 
the estimated residuals of the samples' uncertainty into account (see Fig.[2|. The diffusivity Dz 
and convective velocity RVz are then given by fitting the a cquir ed fiuxes to Eq. (l7| , where after 



the peaking factor is obtained as PF — —^jj^ (see section II.A). 

The instabilities causing the transport are fuelled by the tree energy present in gradients in 
the system, and in general the steeper the gradient the more free energy is available, which is 
expected to lead to stronger modes and more pronounced transport. Two families of gradients 
are available that can drive the instabilities: the temperature gradients {—KVTj/Tj Ki R/Lt) 
and the density gradients {—RVrij/nj w R/Ln.), where j = i, e for main ions and electrons 
respectively]^ Numerical studies have been performed, focused on the dependence of the peaking 
factor on these gradients. 

The main parameters used in the simulations are summarised in Tab. IT] The parameters where 
chosen to represent an arbitrary tokamak geometry at about mid radius, and do not represent any 
one particular experiment. As can be seen in the table, a TE or an ITG mode dominated plasma 
was studied by choosing a steep electron temperature gradient {R/Lxz = 7.0) together with a 
moderate ion temperature gradient {R/ Lt^ = 3.0) to prompt TE mode dominated dynamics. 



'^'the density and temperature gradients of the impurity species can also drive turbulent transport, however, 
for trace amounts this effect is negligible 



and the other way around for ITG mode dominance. It should be noted in this context, that 
TE modes can also be driven by steep density gradients. This case is omitted here, and left 
for future study. In order to preserve quasi-neutrality, Eq. (l6|, Vrig = Vn^ was used. The 
simulations are limited to cases with Tg = Tj. 

In order to ensure that the resolution was sufficient, the resolution was varied separately for 
the perpendicular, parallel and velocity space coordinates, and the effects of this on the mode 
structure, k± spectra and flux levels were investigated. The resolution was then set sufficiently 
high for the effects on the these indicators to have converged. For a typical NL simulation for 
main ions, fully kinetic electrons, and one trace species, a resolution of n^. x n^ x n^ = 96 x 96 x 24 
grid points in real space and of Uy x n^ = 48 x 12 in velocity space was chosen. For QL GENE 
simulations the box size was set to n^: x riy x riz == 8 x 1 x 24 and n„ x n^t = 64 x 12 respectively. 

Simulations have been performed with both Deuterons and protons as main ions, but no 
significant differences in the impurity transport were found between the two cases. 

The impurities were included self-consistently as a third species in the simulations, with the 
trace impurity particle density nz/n^ — 10~^ in order to ensure that they have a negligible effect 
on the turbulence. 

In the present study, a simple s-a geometry is assumed for the simulation domains. The effects 
of different tokamak geomet ries o n drift wave turbulence have been studied in both fluid'^^Mfl 
and gyrokinetic descriptions. '^^'^ 

IV Results and discussion 

For the scalings studied, the charge number Z of the impurities was varied from Z = 2 to Z = 74, 
with a mass to charge ratio AjZ = 2. The scalings of the peaking factor with the temperature 
gradients were studied by varying R/Lt^ between R/Lt^ = 6.0 and R/Lt^ = 10.0 for the 
TE mode case, and similarly by varying R/Lt^z between R/Lt^ z = 6.0 and R/Lt^ z — 10-0 
for the ITG mode case. The density gradient scalings were obtained by varying R/Ln^ between 
R/Ln, = 0.5 and i?/L„, = 5.0. 

QL and NL scalings of PF — —RVz/Dz were obtained using GENE and compared to results 
obtained from the fluid model. 

IV. A Scalings with impurity charge 

The Z scalings of the impurity peaking factor for the TE mode dominated case are presented 



in Fig. 4(a) A good agreement between fluid and NL gyrokinetic results is observed for the 

value kgps — 0.2 used in the QL and fluid simulations. The peaking factors are larger and the 

trends are more pronounced in the QL GENE results, which overestimate the peaking factors 

by approximately a factor of two for low Z impurities. As expected from the discussion in 

section [ll. A| above, PF varies the most for low Z impurities where the thermopinch is stronger. 

For heavier elements, the peaking factor saturates at levels well below neoclassical predictions, 

as seen in previous gyrokinetic and fluid studies, of both TE and ITG mode dominated impurity 
transport.CfiEsIinHlIIini 

For comparison, the results for the ITG mode dominated case is shown in Fig. |4(b)[ The 
two cases show a qualitative difference, with PF falling towards saturation as Z is increased 
for the TE mode case, while the opposite holds for the ITG mode case. This is in line with 
previous QL kinetic and fluid results. ^^"^^1211221 xhe peaking factor is close to zero for low Z 
values in the ITG mode dominated case, however, the sign of PF remains positive for all Z in 
both the TE and the ITG mode dominated case considered. This indicates that a net inward 
pinch is the most common situation in both TE and ITG mode driven impurity transport, for the 



parameters considered. It is, however, known from QL as well as NL gyrokinetic simulations that 
the convection of the impurities can reverse its direction, if the electron heat flux significantly 
exceeds the ion heat fluxJ ^^ I ^^ J 

The qualitative difference between the Z scalings for the TE and ITG mode dominated cases 
can be understood from the balance of the thermodiffusion and parallel impurity compression in 
Eq. ([9]), the two terms having opposite signs for TE and ITG, as discussed above (section II. A I 



The parallel impurity compression is almost independent of Z^ so it can be assumed that the 
thermodiffusion is the main contributor to the observed trends. The thermodiffusion, on the 
other hand, has the strongest effect for low Z values, explaining the drop and rise of PF with 
Z for the TE and the ITG mode respectively. Since this term goes to zero for large values of Z, 
this also explains the observed saturation. 

In the fluid treatment, a strong ballooning eigenfunction is assumed with fcu ^^ = (Sg^i?^) f^ 
Since the contribution from the parallel compression pinch depends on the mode structure along 
the field line, the results are expected to be sensitive to the choice of fcy. To investigate the 
sensitivity of the fluid results to the mode structure, a simplified treatment was used, varying 
ku around its strong ballooning value while keeping the eigenvalues fixed. The results are shown 
in Fig. |5]for kgps = 0.2 and 0.3 in the TE and ITG mode dominated cases. As observed, the 
peaking factors for TE mode turbulence is sensitive to the choice of fcy , with the peaking factor 
going from PF « 2 to PF « when k?, is varied from 0.5 to 2 times its strong ballooning value. 

As is evident from Fig. |4] the value of PF is also dependent on the choice of kgps, the 
perpendicular length scale. Finding the kgps that allows the QL gyrokinetic and fiuid models to 
best capture the behaviour of the impurity transport is non-trivial. For the cases considered, the 
results were obtained with kgps = 0.2. This is in line with previous results regarding comparisons 
of fluid and NL gyrokinetic results.!^ The nonlinear spectra for the fluctuations in the background 
electrostatic potential (</>) are illustrated in Fig. [6] for the TE and ITG mode dominated cases in 
Fig.g 

The spectra both show a peak in the fluctuations at kgps ~ 0.15, well below the wave number 
of maximum linear growth rate, kgps ~ 0.3. We have confirmed that for kgp^ in the range 
0.15-0.4, qualitatively similar QL results are obtained. In the following, kgps = 0.2-0.3 will be 
used. 

A further complication that arises when studying TE mode turbulence is the onset of electron 
temperature gradient (ETG) driven modes. These are mostly sub-dominant, and so are not 
captured by the QL treatment, but may give a nonlinear contribution through the nonlinear 
coupling between the small scale ETG modes and t he lon ger wave lengths of the dominant TE 
modes, and care has to be taken to avoid this effect .'^S'^ 

IV. B Scalings with the temperature gradients 



The obtained scalings of PF with the electron temperature gradient are presented in Fig. 7(a) 
We note that the QL gyrokinetic simulations overestimate the peaking factors by up to ~ 50% 
The fluid results are in good agreement with the NL GENE results. Only weak trends were 



observed, in compliance with previous studies .li^IlTEll As with the Z scaling in Fig. 4(a) the 



NL trend is less pronounced, reaching saturation for lower values of R/Lt^ than the other two 
models. 

For comparison, the results for the ITG mode dominated case are shown in Fig. |7(b)[ As was 
observed for the Z scaling in section [IV.A| above, the trends for the TE and ITG mode dominated 
case are reversed; PF rises with driving gradient for the TE case, but falls for the ITG case. 
The difference between the two trends can be understood in part from the thermodiffusion in 
Eq. ([9|. This term grows more important as the ion/impurity temperature gradient steepens. 



providing a strong outward pinch for the ITG mode dominated impurity transport and thus 
yielding lower values of PF as R/Ltz increases (Fig. |7(b)| . Since the impurity temperature 
gradient is constant for the VTe scaling, however, other effects are behind the TE mode scaling 



in Fig. 7(a) The eigenvalues, in particular the mode growth rates, grow with VTg i, as shown in 
Fig. 7(c) This will alter the relative contributions of the convective terms in Eq. [OJ and hence 



affect the peaking factor. We note here that the eigenvalues in Fig. 7(c) [are normalised to Cg/R 



giving ojr < for TE modes and w^ > for ITG modes. 

As with the Z scaling, the sign of PF usually remains positive for the VTe^^ scalings, though 
a modest flux reversal is observed when the trends of the scalings with Z and R/Lt^ for the 
ITG mode combine. This is the case for He in Fig. |7(b)[ The flux reversal is observed only for 
very steep temperature gradients for the considered parameter values with T^ = Ti. 

IV. C Scalings with density gradient 

In experimentally relevant situations where the impurity and main ion fuelling originates from 
the edge, the core impurity and background density peaking factors should be calculated self- 
consistently for zero particle flux. For this purpose, the equations Pz = and Pi^e = need 
to be solved self-consistently. This is in the following achieved by varying the main ion density 
gradient i?/L„^ until -Tg = is obtained, and using the zero flux background density gradient 
in the impurity transport calculations. We assume trace levels of impurities and use the fluid 
model for simplicity. The results are illustrated in Fig.lSlwhich shows the impurity peaking factor 
R/Ln^ versus R/Ln^ for both the TE and ITG mode dominated cases. The value of R/Ln^ for 
zero background particle flux is marked in the figure. We note that the background density 
peaking is larger than the impurity peaking with R/Ln^ — 3.0 for the TE case and R/Ln^ = 2.5 
for the ITG case. We emphasise that the result is obtained using a collision-less model. It is 
known that collisions have a large impact on the background density peaking in both fluid'221 and 
gyrokinetic models .1^ 

For the i?/L„,, scaling, the same trends are observed in both GENE and fluid data, with a 
strong sensitivity for lower Z impurities. This is particularly evident for the ITG mode case in 
Fig. |8(b)] where the peaking factor for the He impurity shows a marked increase as Vrig steepens 
for both GENE and fluid results, whereas for the heavier elements a nearly flat dependence is 
observed. 

As shown in Fig. 8(c)[ the eigenvalues vary with the electron density gradient. A reduction 



of \ujr\ and an increase of 7 are observed with increasing R/Ln^, which leads to a reduction of 
the relative amplitude of the thermopinch in Eq. (|9| . This explains the observed PF scaling for 
the TE and ITG mode driven cases in Fig. 8 (a) [and Fig. |8(b)| respectively. 



As with the VT^ scaling, the combined effect of the Z and Vue scalings is observed to lead to 
a flux reversal for the He impurity in the ITG mode dominated case in Fig. |8(b)[ This happens 
for flat electron density profiles in the QL GENE results. Outside of this regime the sign of PF 
remains positive. 

V Conclusion and future work 

In the present paper the transport of impurities driven by trapped electron (TE) mode driven 
turbulence has been studied. Non-linear (NL) gyrokinetic simulations using the code GENE were 
compared with results from quasilinear (QL) gyrokinetic simulations and a computationally 
efficient fluid model, viable for use in predictive simulations. The main focus has been on 
model comparisons for electron temperature gradient driven turbulence regarding the sign of the 
convective impurity velocity (pinch) and the impurity peaking factor (R/Ln^) for zero impurity 



flux. In particular, the scaling of the peaking factor with impurity charge Z and with driving 
temperature gradient has been investigated and compared with the more well known results for 
Ion Temperature Gradient (ITG) driven turbulence. 

For the scaling of the peaking factor with the impurity charge Z, a weak dependence was 
obtained from NL GENE simulations, which was reproduced well by the fluid simulations. The 
QL GENE results showed a stronger dependence for low Z impurities and overestimated the 
peaking factor by up to a factor of two in this region. As in the case of ITG dominated turbulence, 
the peaking factors were found to saturate as Z increased, at a level much below neoclassical 
predictions. However, the scaling with Z was found to be weak or reversed as compared to the 
ITG case, where the larger peaking factors were obtained for high Z impurities. 

Using the fluid model it was shown that the impurity peaking factors in the TE mode dominate 
case are sensitive to the mode structure along the fleld lines (fc|| ) through the parallel compression 

pinch. It was shown that assuming a strong ballooning eigenfunction with fc^ = (Sg^i?^) , 
together with fcgps = 0.2, gave a good agreement with the results from the NL GENE simulations. 

The scaling of impurity peaking with the driving background temperature gradients were 
found to be weak in most cases. The QL results were also here found to significantly overestimate 
the peaking factor for low Z values. 

The main ion peaking relative to the impurity peaking was studied using a self-consistent 
treatment of the main ion and impurity particle fluxes. It was found that the main ion peaking 
was slightly larger than the impurity peaking, for both TE and ITG mode dominated turbulence. 
These results were obtained using the fluid model in the collision- less limit. 

The present study is based on low /? plasmas in a simple s-a circular tokamak equilibrium. 
Future work will aim to study the effects of more realistic geometries, flnite /3, as well as effects 
of plasma rotation on impurity transport in NL fluid and gyrokinetic descriptions. 
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TAB. I: Parameters used in the gyrokinetic simulations, ' denotes scan parameters 
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1.0 
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FIG. 1: Impurity flux (/z) dependence on the impurity density gradient (—KVnz/nz = R/Ln^), 
illustrating the peaking factor (PF), the diffusivity (Dz) and pinch (RVz), and the validity of the 
linearity assumption in Eq. (l7| of Fz for trace impurities. Parameters of Eq. (|7| are estimated 
from the calculated fluxes, taking the estimated error of the data into account. The flux is 
acquired as the average of a time series after convergence, as is illustrated in Fig. [21 Data from 
NL GENE simulations of TE mode driven turbulence with He impurities and parameters as in 



Fig. 4(a) The error bars indicate an estimated error of one standard deviation. 
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FIG. 2: Time series showing fluctuations in the main ion density (nj^) and impurity flux (r^) after 
averaging over the whole flux tube (see Fig. Isl. The averaged impurity flux {{Fz)) is calculated 
from Fz, discarding the flrst portion to ensure that the linear phase of the simulation is not 
included. {Fz) is used for flnding the peaking factor for the impurity species, as is illustrated in 
Fig. IT] The bursty nature of the transport is seen in the peak around t « 185 R/cg. These bursts 
have been found to affect the average flux little, but to signiflcantly increase the estimated error in 
{Fz) ( ). Data from NL GENE simulation of TE mode driven turbulence with He impurities. 



The parameters are the same as in Fig. 4(a) with —RVnz/nz = R/Lm = 1-5 
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FIG. 3: A cross-section of the flux tube, showing the fluctuation of the electrostatic potential 
(f>. Data from NL GENE simulation of TE mode turbulence, with parameters as in Fig. 4(a) at 
t« 300 i?/c,.. 
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4(a): dependence of the peaking factor (PF) on Z for the TE case 
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4(b): dependence of the peaking factor {PF) on Z for the ITG case 

FIG. 4: Scalings of the peaking factor (PF) with impurity charge (Z). Parameters are q = 1.4, 
s = 0.8, e = r/R ^ 0.143 in b oth su bfigures, with R/Lt, = R/Lt^ = 3.0, R/Lt, = 7.0, 
i?/L„^ = 2.0 for the T E case (Fig. [4(a)| ) , and R/Lt, = R/Lt^ = 7.0, R/Lt, = 3.0, R /L^^ = 3.0 
for the ITG case (Fig. |4(b)| ). The error bars for the NL GENE results in Fig. 4(a) indicate an 
estimated error of one standard deviation. 
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FIG. 5: Scaling of the peaking factor (PF) with kf./kf. ^^^ for He impurity, where k 



3q R ) is the strong ballooning value; fluid results with parameters as in Fig. 4(a) (TE) 



and Fig. 4(b) (ITG) 
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Fig. 4(a) (TE) and Fig. 4(b) (ITG) 
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7(a): dependence of the peaking factor (PF) on the norniahsed electron temperature gra- 
dient for the TE case 
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7(b): dependence of the peaking factor {PF) on the normalised ion temperature gradient 
for the ITG case; 
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7(c): real frequency (cj,.) and growth rate (7) for the two cases in Fig. 7(a) and Fig. 7(b) 



FIG. 7: Scalings of the peaking factor (PF) with the electron and ion temperature gradients 
{-RVTe^i/T^^i = R/Lt,J- Para meter s for the TE and ITG mode case as in Fig. [i] with 
kePs = 0.2. The eigenvalues in Fig. 7(c) are from QL GENE simulations, they are normalised to 

Cs/R. 
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8(a): dependence of the peaking factor {PF) on the nomiaUsed electron density gradient 
for the TE case, also indicated is the main ion peaking factor (PF^) from fluid theory 
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8(b): dependence of the peaking factor {PF) on the normalised electron density gradient 
for the ITG case, also indicated is the main ion peaking factor (PFe) from fluid theory 
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8(c): real frequency (cj,.) and growth rate (7) for the two cases in Fig. 8(a) and Fig. 8(b) 



Scalings of the peaking factor {PF) with the electron density gradient {—RVne/rie = 
Parameters for the TE and ITG mode cases as in Fig. [4] with kgpg — 0.3 for both cases. 



FIG. 

The eigenvalues in Fig. 8(c) are from QL GENE simulations, they are normalised to Cg/R. 
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